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Abstract 

The influence of lateral adsorbate diffusion on the dynamics of the first-order phase 
transition in a two-dimensional Ising lattice gas with attractive nearest-neighbor 
interactions is investigated by means of kinetic Monte Carlo simulations. For ex- 
ample, electrochemical underpotential deposition proceeds by this mechanism. One 
major difference from adsorption in vacuum surface science is that under control 
of the electrode potential and in the absence of mass-transport limitations, local 
adsorption equilibrium is approximately established. We analyze our results using 
the theory of Kolmogorov, Johnson and Mehl, and Avrami (KJMA), which we ex- 
tend to an exponentially decaying nucleation rate. Such a decay may occur due to 
a suppression of nucleation around existing clusters in the presence of lateral ad- 
sorbate diffusion. Correlation functions prove the existence of such exclusion zones. 
By comparison with microscopic results for the nucleation rate I and the interface 
velocity of the growing clusters v, we can show that the KJMA theory yields the 
correct order of magnitude for Iv 2 . This is true even though the spatial correlations 
mediated by diffusion are neglected. The decaying nucleation rate causes a gradual 
crossover from continuous to instantaneous nucleation, which is complete when the 
decay of the nucleation rate is very fast on the time scale of the phase transfor- 
mation. Hence, instantaneous nucleation can be homogeneous, producing negative 
minima in the two-point correlation functions. We also present in this paper an re- 
fold way Monte Carlo algorithm for a square lattice gas with adsorption/desorption 
and lateral diffusion. 
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1 Introduction 

The electrochemical adsorption of a two-dimensional film on a metal elec- 
trode can involve a phase transition at the surface [1]. When the adsorbate- 
adsorbate interactions are attractive, this phase transition is of first order; 
hence, below a critical temperature, there is a discontinuity in the coverage at 
a certain electrochemical potential. Stepping across the electrostatic potential 
difference between electrode and solution that corresponds to this coexistence 
point induces adsorption (or desorption). Whereas, thermodynamically, this 
phase transformation is abrupt rather than gradual, its kinetics are governed 
by the microscopic mechanism by which it proceeds. On a flat, defect-free 
single-crystal surface, the decay of the metastable phase takes place through 
nucleation and growth of two-dimensional clusters. The classical theory for 
this kind of phase change was developed in the 1930s and '40s by Kolmogorov, 
Johnson and Mehl, and Avrami (KJMA theory) [2-6]. Finding successful ap- 
plication in many fields of science, it is also the basis for the interpretation 
of kinetic data of electrodeposition processes that proceed by nucleation and 
growth [7,8]. However, the assumptions of the theory may not always be 
met in real experimental systems. In particular, an earlier study from our 
group indicates that homogeneous nucleation can also be instantaneous, and 
that there might be a gradual crossover from continuous to instantaneous 
nucleation [9]. This casts some doubt on whether the nucleation is spatially 
uncorrelated [10,11]. 

Lattice-gas models are now a widely applied tool in the study of electro- 
chemical adsorption [12,13]. They are appropriate when the adsorbate layer 
is commensurate with the substrate [14]. When the adsorbate is in quasi- 
equilibrium with the bulk under control of the electrode potential, a grand- 
canonical Hamiltonian should be used. Kinetic Monte Carlo simulations [15] 
of lattice-gas models can be used to study the dynamics of such systems when 
they can be well approximated as a superposition of independent, stochastic 
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processes [12]. This is generally true when the elementary processes evolve 
through thermally activated crossings of well-defined transition states. Then, 
one should choose an Arrhenius-type Monte Carlo dynamic [16, 17]. For the 
investigation of the model dynamics at macroscopic time scales this is the 
method of choice, since it does not rely on mathematical approximations and 
since it coarse-grains over the fluctuations (vibrations) whose time scales are 
orders of magnitude faster than adsorption, diffusion, etc. With this method 
real problems have been tackled and successfully described [18], for example 
for underpotential deposition of metals [19,20]. However, the problem may 
require substantial effort in providing a sufficiently realistic description. For 
example, one may need to include coadsorption of anions, long-range or many- 
body potentials [14,18,21], etc. 

Here we follow a somewhat different approach. Rather than in a specific sys- 
tem, we are interested in the general behavior of two-dimensional first-order 
phase transitions. Hence we choose as a model system a two-dimensional 
square lattice gas with nearest-neighbor attractive interactions. No solvent 
or counter ions are included. This model can be interpreted as a simplified 
description of electrochemical metal underpotential deposition or of the early 
stages of epitaxial growth when no nucleation takes place in the second layer. 
However, since it is isomorphic to the square-lattice Ising model [22], it is far 
more general and may, among other phenomena, also describe ferromagnetic 
systems [23]. The focus of this paper is on the influence of nearest-neighbor 
adsorbate diffusion on the dynamics of the phase transformation. Unlike in vac- 
uum surface science, adsorbate diffusion is not a condition for nucleation since 
fluctuations in cluster sizes can also be mediated by the quasi-equilibrium with 
the solution (instead of the two-dimensional gas of adsorbed particles). Like- 
wise, an earlier study [23] has shown that the phase-transformation dynamics 
of this lattice gas can be well described by the K JMA theory for continuous nu- 
cleation over a wide range of potentials in the absence of adsorbate diffusion 
(in magnetic terms, spin flip only). Another preliminary study [9] suggests 
that in the presence of diffusion (spin exchange in magnetic language) the 
dynamics may change. For fast diffusion, the phase transformation appeared 
to proceed by instantaneous nucleation, and for intermediate diffusion rates 
there was a gradual crossover. This was interpreted as a result of exclusion 
zones around the clusters, in which nucleation was suppressed. It was specu- 
lated that these exclusion zones may become space-filling in the early stages 
of the phase transformation. The crossover also changed the morphology of 
the adsorbate phase. However, the study used Glauber instead of Arrhenius 
dynamics and was limited to one value of the electrochemical potential. Here 
we extend the study to a wide range of potentials. Furthermore, we present a 
more stringent analysis of the kinetic data and use microscopic details of the 
simulation for comparison with the theory. 

The remainder of this paper is organized as follows. In Sec. 2 we review the 
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KJMA theory and extend it to an exponentially decaying nucleation rate. In 
Sec. 3 we introduce the lattice-gas model and the details of the algorithm and 
the model calculations. Results are presented thereafter. In Sec. 4.1 we describe 
the phase transformation dynamics and the morphology of the adsorbate phase 
and give kinetic parameters according to the theory. In Sec. 4.2 we present 
the results for the microscopic rates. In Sec. 4.3 we discuss the accuracy of the 
continuum description of the KJMA theory, and in Sec. 4.4 we discuss under 
which conditions the crossover to instantaneous nucleation occurs. A summary 
and conclusions are given in Sec. 5. A description of the n-fold way algorithm 
for adsorption/desorption and diffusion can be found in the Appendix. 



2 Theory 



We here briefly review the KJMA theory and describe how we adapt it for a 
gradual crossover between continuous and instantaneous nucleation. 

The basic assumption of the KJMA theory is that nucleation events are ran- 
dom and uncorrelated. Then, it is mathematically quite simple to correct for 
the overlap of growing circular droplets, provided that they do not change 
shape upon coalescence [24,25]. One defines an extended coverage 9 cxt by 
summing up the areas of all growing droplets regardless of overlap with other 
clusters, and normalizes this by the total surface area. When the radius r of 
a cluster is a function of its birth time t, this reads 

/•fmax of* ~ o ~ 

#ext = / p(r)irr 2 dr = nv 2 / I(t)(t - t) 2 dt : (1) 
Jo Jo 

where p(r) is the number of clusters of radius r per surface area, and I(t) is 
the nucleation rate at time t. The last equality is for circular clusters with a 
radius growing linearly with an interface velocity v. It is important to realize 
that this nucleation rate includes nucleation events in the transformed area, 
since otherwise the condition of random, uncorrelated nucleation would be vi- 
olated [26,27]. There has been some confusion concerning this point [28]. We 
therefore call the nucleation rate I(t) the extended nucleation rate wherever 
the distinction is important. In the system we consider in our investigation, 
nucleation is homogeneous and proceeds by thermal fluctuations. It can be de- 
scribed by classical nucleation theory [22], but we do not address this question 
further here. 

When I is constant throughout the phase transformation, we speak of contin- 
uous nucleation, and performing the integration in Eq. (1) yields 

nlv 2 t 3 

#cxt = • (2) 
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Here, / has the dimension l/(areaxtime). It is clear that the "real" nucleation 
rate I T is reduced by the fraction of the covered surface 9 and eventually de- 
creases to zero; this is not to be confused with a decaying extended nucleation 
rate (see below). At early times, though, the real nucleation rate can be used 
as an estimate for the extended nucleation rate, since 6 ps 0. 

In contrast, when all nuclei are already present at the very beginning of the 
phase transformation and no further nucleation takes place, I(t) = IS(t), 
where S(t) is a Dirac delta function, and the integration in Eq. (1) yields 

9 ext = nlvh 2 . (3) 

This instantaneous nucleation is often associated with heterogeneous nucle- 
ation at defects. Here, / has the dimension 1/area since it is virtually a surface 
density of nucleation centers. 

The real coverage can be calculated from the extended coverage 9 by [5] 

l-^ = 0KJMA = e^ xt . (4) 

The quantity 0kjma is called the relaxation function (in the KJMA theory). 
We will later define a relaxation function for our model system and use 0kjma 
as a coarse-grained approximation for it [23]. 

2. 1 Extended volume for decaying nucleation rate 

The gradual crossover from continuous to instantaneous nucleation with in- 
creasing diffusion rate implies that, for some intermediate diffusion rates, the 
extended nucleation rate / is neither a constant nor a delta function, but de- 
cays with time. In the following, we deduce how an exponentially decaying 
extended nucleation rate 

I(t) = Ioe- Xt (5) 

with Iq the initial nucleation rate and A an inverse decay time translates into 
a 9{t) relation for the coverage. The dimension of I is l/(area x time). 

In Eq. (5), I(t) is independent of the position. This is a mean-field approxima- 
tion, in which any local reduction of the nucleation rate is smeared out over 
the whole surface. Within this approximation, nucleation remains spatially 
uncorrelated, and Avrami's law [Eq. (4)] holds. Hence we need to calculate 
the extended coverage that follows from Eq. (5). Inserting I(t) into Eq. (1) 
and performing the integration yields 

«-W = ^ - ^ + ^ (1 - e-») . (6) 
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In the limit of At ^> 1, the second and the third term are of the order of At 
and A 2 t 2 smaller than the first term, and we recover the well-known result for 
instantaneous nucleation, Eq. (3), with I = In the limit At <C 1, the ex- 
ponential can be expanded up to third order, which is the lowest non-vanishing 
order, and we recover the well-known result for continuous nucleation, Eq. (2), 
with I = I . 

We note that the exponential decay of the extended nucleation rate is not de- 
rived from a physical model, but rather simply a mathematically convenient 
assumption. Moreover, in real systems as well as in our microscopic model 
system, the reduction of the nucleation rate is thought to be local in exclusion 
zones around existing clusters. This would violate the assumption that nu- 
cleation events are uncorrelated. Such correlations are ignored by the theory 
presented. A theory that includes these effects could not make use of Avrami's 
law. 

A fit of experimental or simulation results for the relaxation function to 
Eqs. (4) and (6) can yield estimates for A and for Iov 2 , but no separate esti- 
mates for J and v. 



3 Model and algorithm 

3.1 Model 

We model the electrodeposition process using a square lattice-gas model with 
attractive nearest-neighbor interactions. Since the subcritical fluctuations of 
the adsorbate are at quasi-equilibrium, a grand-canonical Hamiltonian is used: 

W = -S -JiEci(f), (7) 

(i,3) 

where $ is the lateral interaction constant and \x the electrochemical potential 
of the adsorbing ion in the solution 2 . Mass-transfer limitations from the 
solution to the surface are neglected. The occupation variables q take the 
value when site i is empty, and 1 when it is occupied. 

We perform simulations for an L x L = 256 x 256 lattice with periodic bound- 
ary conditions, using for all simulations the interaction constant <£> = 4 and the 

2 For a spin s« = ±1 Ising model with z nearest neighbors (here, z = 4) and J\f 
sites, coupling constant J and external magnetic field H, the mapping between the 
models is Hi = H LG + (Af/2) [n - (/x /2)], a = ( 8i + l)/2, $ = 4 J, n - no = 2H, 
and = ~2zJ [22]. 
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temperature T = 0.8 T c pa 0.454$ = 1.815, where T c is the critical tempera- 
ture of the Ising model for the given parameters [29] . Energy and temperature 
units are chosen such that Boltzmann's constant k B = 1. 



3.2 Algorithm 



The dynamics of the phase transformation are investigated using kinetic Monte 
Carlo simulations [15]. We assume that the model dynamics can be reason- 
ably approximated by a few stochastic elementary processes; hence, we choose 
one-step Arrhenius-type dynamics [16, 17] for the transition probabilities, in- 
troducing a local transition barrier A. The transition rate R for a transition 
from configuration a to configuration b is approximated by the corresponding 
rate for a transition from a to a fictitious transition state: 

R(a ^b) = u exp I — J , (8) 

where TO is the energy at the transition state. Using a symmetric (a = 1/2) 
Butler- Volmer- (Br0nsted-)type expression for TO [12]: 



A \ / H(b)-H(a) \ 
R(a ^b) = u exp -— - exp — — , (10) 



one obtains 



k B TJ r \ 2k B T 

where z/ sets the Monte Carlo time scale and will be chosen as unity in the 
following. Another choice for an Arrhenius-type dynamic is the transition- 
dynamic algorithm [30,31]. 

For each kind of elementary step that is included in our model simulations, we 
specify a value for the transition barrier A. In particular, these are A ad / dcs for 
adsorption/desorption steps and A dif for diffusion steps to a nearest-neighbor 
site. The former is equivalent to a spin-flip step, the latter to a spin-exchange 
step in the terminology of magnetic spin models. By varying the difference 
A ad/dcs — A d if we can set the ratio between adsorption/ desorption and diffusion 
steps. The absolute values have no particular meaning unless we attempt to 
specify the relation between physical and Monte Carlo time scales [18]. We 
treat A as merely formal and note that, according to Eq. (9), TO might not 
be larger than both Tt{a) and Tt{b) for a strongly exothermic or endothermic 
step, viz not represent a true barrier, even for positive A. In the Monte Carlo 
algorithm chosen, transition rates R > 1 pose no problem (see below). 
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At high diffusion rates, standard Monte Carlo algorithms are prohibitively 
slow since much simulation time is spent on the coverage- conserving diffu- 
sion steps. It is therefore necessary to use more advanced methods: we choose 
the n-fold way algorithm of Bortz, Kalos and Lebowitz [32] and Gilmer [33]. 
This algorithm is rejection- free: of all possible new configurations, one is cho- 
sen with probability Rj/R to t, where R to t = an d then accepted with 
probability 1. No time is spent on choosing sites and calculating energy dif- 
ferences for rejected moves. After each step the simulation time is updated 
by — (l/i?tot) InrMCSS (Monte Carlo steps per site), where r is drawn from 
a uniform pseudo-random distribution between and 1. We need to keep a 
list of the transition rates for all possible moves from the actual configuration. 
This is done for a single-site n-fold way as follows. For each site s in the lat- 
tice we compute the transition rate R s , m f° r every possible elementary step m 
according to Eq. (10): adsorption/desorption (R St u adsorption when site is 
empty, desorption otherwise); exchange with nearest neighbor to the right 
{R s ,2\ equal to zero when both are occupied or both are empty), exchange 
with nearest neighbor below (R s ,3', equal to zero when both are occupied or 
both are empty). This minimizes the possible steps per site and makes sure 
that every possible diffusion step on the lattice is counted once. We then cal- 
culate the total rate at which site s performs a move as S s = J2 m Rs,m- I n a 
Monte Carlo step, a random number is drawn from a uniform distribution to 
choose with probability S s /R to t the site s that performs the next move. Then, 
another random number is drawn to choose with probability R SjTn / S s which 
elementary step m is performed. Finally, a third random number is drawn to 
choose the time step by which the simulation clock is updated. After that, the 
respective R and S of all sites that are affected by the move are updated. In 
order to accelerate this process, the latter are stored in a binary tree. 



Recognizing that there is a finite number of configurations of a site plus en- 
vironment, each characterized by an energy that is constant throughout the 
simulation, it is in principle possible to tabulate the energy differences that are 
connected with each change of local configuration (called classes). This has pre- 
viously been worked out for the square Ising lattice without spin-exchange (dif- 
fusion) [32], where 10 classes are to be considered for adsorption/desorption. 
Here, we additionally have to define 32 classes for the diffusion step; our n-fold 
way is hence a 42-fold way algorithm. The situation is further complicated by 
the fact that the diffusion step calls for including a larger portion of the lat- 
tice into the affected environment of a site. We have to consider a total of 21 
classes up to third-nearest neighbor positions after an adsorption/desorption 
step and of 31 classes up to fourth-nearest neighbor positions after a diffusion 
step. Details are shown in the Appendix. 
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3.3 Model calculations 



Our goal is the investigation of the influence of diffusion on the dynamics of 
the first-order phase transition. To that end, we vary the diffusion barrier A^f 
between 5 and 15, keeping A a( j/des — 15 constant. Moreover, we perform simu- 
lations without diffusion, corresponding to Adif = oo. We start our simulations 
with an empty lattice, and at t = we set the electrochemical potential to a 
value between \i — /j, = 0.4 and /i — /i — 1.6, where /j,q is the electrochemical 
potential at coexistence. For all these values of the electrochemical potential, 
the phase transformation is in the multi-droplet regime [23], where the num- 
ber of supercritical adsorbate clusters is large. For each set of parameters, we 
perform 500 to 1000 independent simulation runs. This high number is neces- 
sary to accurately obtain the variance of the coverage. For quantities that do 
not depend on higher moments, we typically take statistics from only about 
100 runs. 

We sample the coverage 9 as a function of the simulation time. Moreover, we 
follow the size of every cluster in the lattice throughout the simulation until 
the coverage reaches = 1/2; this enables us to get microscopic interface 
velocities and nucleation rates. Later, cluster sizes are obtained using the 
Hoshen-Kopelman algorithm [34] at certain coverages. We compute the time- 
dependent two-point correlation function (see Ref. [23], correcting for missing 
factors of powers of L) 

G(r,t) = ( Ci (t) Cj (t))-m)) 2 , (11) 

where r = f— fj, as 1/ ^/T? times the inverse Fourier transform of the structure 
factor 

S(q,t) = (c q ^)-L 2 (9(t)) 2 5^, (12) 
where 5^$ is the Kronecker delta function and is the Fourier transform of 

Cf 

1 L2 

Cq(t) = -j= Sci(t) exp(-i(f • fl). (13) 

i= l 

The angular brackets denote averages over independent runs. Due to the use 
of a Fast Fourier Transform algorithm [35], this approach is significantly faster 
than the direct calculation of G(r,t). We only report the circularly averaged 
G(r,t). 



4 Results 

Our interest is in the influence of diffusion on the kinetics and dynamics of 
the phase transformation. It is therefore necessary to first characterize how 
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the fraction of steps that are diffusion steps depends on the model param- 
eters. We count the diffusion steps and divide their number by the number 
of total steps. Since the time step by which the clock is updated after a suc- 
cessful move does not depend systematically on the type of move performed, 
this simpler approach is approximately equivalent to computing true rates of 
the elementary processes. Figure 1 shows that the fraction of diffusion steps 
changes only weakly during the phase transformation. Hence, it is sufficient 
to discuss the average fraction of diffusion steps throughout the whole simu- 
lation. In Table 1, we show that the diffusion fraction increases strongly with 
decreasing diffusion barrier A dif , reflecting the exponential law, Eq. (10), and 
increases weakly with the electrochemical potential, reflecting that diffusing 
particles spend on average longer time on the surface before they desorb. With 
^dif = 5, diffusion is the predominant process. 

The critical cluster size n* is governed by the balance of volume and surface 
free-energy terms of the clusters [22,36]; it hence depends on the interaction 
constant $ (which we keep constant), and on the electrochemical potential /z 
(which we vary). The critical cluster size is the number of particles n in the 
cluster for which the free energy F{n) has a maximum. We calculate F(n) = 
— k-QT\a.Z{n) from the restricted partition function Z(n), which we obtain 
from enumerating lattice animals, hence counting all possible configurations 
of a given cluster size [9]. The resulting critical cluster size is shown as a 
function of \i in Table 2. As expected, it decreases with increasing \x. 



4-1 Dynamics of the phase transformation 



We characterize the phase transformation using the relaxation function <fi [23], 
which is closely related to 1 — 6 and decays from 1 to during the adsorption. 
It is defined as 

«')=?H « i4 » 

where 6 S is the stable coverage, which depends on the electrochemical potential 
and is easily obtained from equilibrium Monte Carlo simulations. We define 
the metastable lifetime r as the time when = 1/2. 

The influence of the electrochemical potential on the kinetics of the phase 
transformation is shown in Fig. 2. The phase transformation lags behind the 
potential switch, but it is markedly accelerated when the electrochemical po- 
tential, which is the driving force for adsorption, is higher. In Fig. 3 we show 
the influence of diffusion on the kinetics for \i — /i = 0.4. At A dif = 15, the 
phase transformation is only slightly faster than in the absence of diffusion; 
however, at lower values of the diffusion barrier, the acceleration is consider- 
able. This trend is verified in Table 3 for all electrochemical potentials and 
diffusion barriers considered. 
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In a previous investigation [9], using Glauber dynamics, we found that when 
diffusion is sufficiently fast on the time scale of adsorption/desorption, there is 
a crossover from continuous to instantaneous nucleation. We demonstrate in 
Fig. 4 that this behavior persists when introducing local barriers into the dy- 
namics. We normalize the time scale with the metastable lifetime; then, when 
only the time scale of the phase transformation changes and not the dynamics, 
the plots should collapse. This is the case for the simulations without diffusion 
and with A dif = 15, but for lower diffusion barriers we find increasingly strong 
deviations. Figure 4(a) shows the plot for the model of continuous nucleation 
[Eqs. (2) and (4)]. The deviations from linearity at early times probably come 
from an overlap of the time scales of the formation of the metastable phase 
(starting with a completely free surface) and of the metastable decay. However, 
in the absence of diffusion, there is a marked region where linear behavior is 
found. This region becomes shorter and eventually disappears when diffusion 
becomes predominant. Instead, as shown in Fig. 4(b), for A dif = 5 the predic- 
tion for instantaneous nucleation [Eqs. (3) and (4)] is fulfilled. Similar results 
were reported in Ref. [9]. 



4-1.1 Crossover to instantaneous nucleation 

On the basis of the results from the previous section, we attempt to fit the 
phase transformation curves to our model for decaying nucleation rates (Sec. 2.1) 
According to Ref. [23] , when the K JMA theory is used as a coarse-grained ap- 
proximation for the microscopic Ising (or lattice-gas) model, the relaxation 
function can be corrected for the fact that the simulation does not start from 
the metastable phase: 

^(*) ~ Tife — tt^kjmaW- (15) 



In a first step, we fit (f)(t) from the simulations without diffusion to Eqs. (2), 
(4), and (15) and obtain estimates for the two free parameters of the fit, 
Iv 2 and 9 ms . To do so, we need an estimate for 9 S , which we obtain from 
equilibrium Metropolis Monte Carlo simulations of the order of 10 6 MCSS. The 
thermodynamic quantity 6 ms is only a function of T and fi and independent of 
the microscopic details of the dynamics. It can therefore be used later for the 
simulations with diffusion at the same potential. Strictly, due to the imperfect 
separation of the timescales for the metastable decay and the formation of the 
metastable phase from the empty surface, the estimates for 9 ms obtained in this 
way are compromised by a systematic error. In particular at higher potentials 
they appear to be too high, and we refrain from reporting their numerical 
values. However, the error is most likely dominated by the potential and not 
by the diffusion rate, and we use the same 9 ms for all diffusion rates at the 
same potential in the following. 
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In a second step we fit <f>(t) from all the simulations (with and without diffu- 
sion) to Eqs. (4), (6), and (15). The free parameters of the fit are I^v 2 and A. 
The former contains information about the kinetics of the phase transforma- 
tion, the latter about the decay time of the nucleation, viz the dynamics of 
the phase transformation. We fit the relaxation function in the interval from 
0.5r to r, since at short times the decay to the metastable phase affects <p, 
and at longer times coalescence effects become noticeable [9]. By dividing the 
data into about 10 to 20 subsets and performing the fit independently for each 
subset, we obtain the standard deviation of the fitting parameters, and hence 
an estimate of the standard error. 

In Table 4 we show the results for the kinetic parameter IqV 2 . The standard 
error is insignificant for the number of digits given. The values of Iv 2 reflect 
the same trends as the metastable lifetime r. In the absence of diffusion and 
for A dif = 15, the values are identical to those obtained for Iv 2 from fits to 
Avrami's law for continuous nucleation. For the values for which we expect 
a possible crossover to instantaneous nucleation to be most complete (// — 
/j,q = 1.6, A^f = 5) Iov 2 /\ = 1.3 x 10~ 9 , in satisfactory agreement with 
Iv 2 = 9.5 x 10~ 10 from a fit to Avrami's law for instantaneous nucleation. 

This behavior is reflected in the values of Ar that we report in Table 5. When 
Ar = 1, the nucleation rate has decayed by a factor 1/e at the metastable 
lifetime. Hence, when Ar 1, we are in a transition regime where the Avrami 
plots neither resemble instantaneous nor continuous nucleation. In the regime 
of continuous nucleation, Ar <C 1; the exact numerical value is insignificant 
and thus scatters. This can cause quite large relative standard errors. When 
Ar becomes larger, we are approaching the regime of instantaneous nucleation. 
These finite values of Ar are indeed significant, as shown by the quite small sta- 
tistical errors in these cases. As the results in the previous paragraph suggest, 
for the highest potential and lowest diffusion barrier considered, the nucleation 
is virtually instantaneous. A value of Ar « 4 means that the extended nucle- 
ation rate I is reduced to about Jo/50 when the relaxation function reaches 
1/2. We show the fits for ji — /i = 1.2 in Fig. 5 to demonstrate how the 
model captures the crossover from continuous to instantaneous nucleation. 
Apart from deviations at early times, the fits are satisfactory. The decrease in 
curvature with increasing diffusion rates is well captured. 



4-1.2 Morphological changes 

In Ref. [9] we demonstrated that the crossover from continuous to instan- 
taneous nucleation is accompanied by changes in the adsorbate morphology 
during the phase transformation. In particular, the cluster-size distributions 
showed a depletion of mid-sized clusters and an enrichment of large clusters. 
Here we show that similar behavior is found using dynamics with local barri- 
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ers. 



According to Sekimoto's theory for uncorrelated, continuous nucleation [37, 
38], the two-point correlation function G(r) is a monotonically decreasing func- 
tion of the distance r. We show the results for /i — /i — 0.4 and A^f = oo in 
Fig. 6 (a). The correlation functions indeed decay monotonically for all cover- 
ages. As can be easily shown from the definition of G(r), G(0) = (9) — (9) 2 = 
{9) (1 — (9)), which equals the mean-square fluctuation of the occupation vari- 
able from its average [39]. For later use, it is convenient to normalize the corre- 
lation function as G(r) = G(r)/G(0), such that (5(0) = 1. The un-normalized 
G(0) assumes the same value for coverages that add up to 1, and is maximal 
for 9 = 0.5. The high-coverage correlation function of these pairs always de- 
cays more slowly, showing the asymmetry of the morphology about 9 = 0.5. 
It has been shown in Ref. [23] that the agreement of Sekimoto's theory with 
G(r) is very good in the absence of diffusion. Only at short length scales are 
there some deviations due to the in-phase correlations which are neglected in 
the theory. The effect of diffusion on G(r) is shown in Fig. 6 (b). Here, for 
/i — no — 1.6 and A dif = 5, the correlation function decays to negative values, 
and after going through a minimum converges to 0. The negative values of 
G(r) reflect the exclusion zones around existing clusters, in which nucleation 
is diminished since aggregation to the cluster is more likely than the formation 
of a new cluster. This result is very important since it shows a violation of 
the assumptions behind the KJMA theory and makes the use of Avrami's law 
in the interpretation of the phase-transformation kinetics questionable. Any 
quantitative errors of the nucleation rate and interface velocity estimated on 
the basis of the KJMA theory are most likely due to these diffusion-induced 
correlations. We note that similar minima are also found in the two-point 
correlation functions during phase separation with a conserved order parame- 
ter [40] , but not during phase ordering with a non-conserved order parameter. 
As in the case studied here, negative correlations during phase separation are 
mediated by diffusion. In our simulations, the minima occur first when At is of 
the order of 1 (see Section 4.1.1). They are most pronounced around 9 = 0.5. 
With increasing coverage they move toward larger r since the average cluster 
size increases. 

From the inverse of the initial slope of the normalized correlation function 
G(r) we can estimate the correlation length /. As we show in Fig. 7, the 
correlation length as a function of the coverage is a concave function with the 
maximum shifted to values larger than 9 = 0.5. This shows that correlations 
are of longer range for the vacancy clusters at coverage 1 — 9 than for the 
corresponding adatom clusters at coverage 9. Generally, the correlation length 
decreases with increasing /i, and increases with decreasing Adif, indicating 
that diffusion causes the formation of larger structures on average. This is 
corroborated by an investigation of the cluster-size distributions. Figure 8 
shows as a typical example the distributions for // — // = 1-6 at 9 = 0.3. 
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With rapid diffusion, the densities for the largest clusters are enhanced at the 
expense of mid-sized clusters. Similar results were found in Ref. [9]. 

Some of these observations can easily be verified from inspecting snapshots 
of the real-space morphology, like the ones shown in Fig. 9. The asymmetry 
of the morphology about 9 = 0.5 is clearly visible comparing the snapshots 
at 9 — 0.3 and 6 = 0.7 for all potentials and diffusion barriers shown. There 
are fewer, but larger vacancy islands at the higher coverage than adsorbate 
islands at the lower coverage, and their shape is markedly more different from 
circular. The average island size at the same coverage increases with decreasing 
potential and with increasing diffusion rate, causing the observed trends for 
the correlation length. 



4-1.3 Separate nucleation rate and interface velocity 

Within the KJMA theory for spatially random, continuous nucleation, Seki- 
moto [37, 38] has derived an expression for the two-point correlation function 
r(r) as a function of the nucleation rate, the interface velocity, and time. Y(r) 
can be used as a coarse-grained approximation for the correlation function 
G(r) from our simulations. From Eqs. (12) and (13) it can be shown that 

S(0)=L 2 v a i(9)=Y / G(fl). (16) 

i 

The last equality is the inverse Fourier transform for q — 0. Hence, fitting 
the simulation results for L 2 var(#) to Sekimoto's expression for T can yield 
separate estimates for the nucleation rate and the interface velocity / and 
v, respectively [23], provided the assumptions of Sekimoto's derivation hold. 
Taking into account the coarse graining and local fluctuations, one obtains [23] 



L 2 var(£) » [9 ms - 9 S ] 2 8™ 2 t 2 0KjMA * [© (^ 3 ) " \ 



2J 

+ ^KJMA&B^Xms + [1 ~ 0KJMa] ^B^Xs- (17) 



The function 



G(x) = f'ye^^dy (18) 
Jo 

is obtained by numerical integration using 



*(v) = I 



arccos(y) — 2yy / l — y 2 + y 3 In [ — — — 



(19) 



Here, Xs and Xms describe the local fluctuations of the stable and the metastable 
phase, respectively. The former is equivalent to the susceptibility in magnetic 
language and can be obtained from equilibrium Monte Carlo simulations [41]. 
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Iv 2 in the argument of is known from the fits to Avrami's law for decaying 
nucleation rate (Section 4.1.1). The free parameters of the fit are the interface 
velocity v and Xms- Knowledge of v and Iv 2 enables the calculation of /; we 
do not report the results for Xms- 

The model fits the data very well for low supersaturation and high diffusion 
barriers. The quality of the fits decreases a little with increasing supersatu- 
ration. For Adif = 5, though, Sekimoto's model is not capable of describing 
the data. Obviously, diffusion causes a violation of its underlying assump- 
tions. In particular, nucleation is no longer spatially random. We therefore 
exclude all simulations which show minima in the two-point correlation func- 
tion G(r) from this analysis. Results for the nucleation rate and the interface 
velocity are shown in Tables 6 and 7, respectively. The nucleation rate grows 
strongly with the supersaturation — by more than two orders of magnitude 
from \i — fi = 0.4 to fj, — Hq — 1.6. The interface velocity increases in the same 
sense, but less dramatically, by less than an order of magnitude. Due to the 
limitations of the model, no statement can be made about the trend with (no- 
ticeable) diffusion rates. Hence, and also for comparison with the data without 
diffusion, we will in the following use information from the microscopic details 
of the simulations. 

4-2 Microscopic rates 

To obtain microscopic nucleation rates and interface velocities, we follow the 
time development of the size (number of particles) n of every cluster during a 
simulation run. An adsorption step can cause the coalescence of two or more 
clusters. A desorption step can cause the dissociation of a cluster; one needs to 
check every time if the possible fragments are still connected. The algorithm 
can be accelerated by realizing that the maximum number of fragments is the 
number of nearest neighbors of the desorbed particle. Moreover, it is faster to 
look for connections in the vicinity of the desorbed particle first. 

We count coalescence and dissociation events. The evaluation of nucleation 
events is described in Section 4.2.2. Diffusion is treated as an adsorption step 
followed by a desorption step. Beyond the percolation threshold the algorithm 
is prohibitively slow, and we use it only up to 6 = 0.5. 

4-2.1 Interface velocity 

The KJMA theory in the present form makes use of the common assump- 
tion that the radius of a circular cluster grows linearly with time (Allen-Cahn 
approximation [22,37,38]). This assumption is expected to be very good for 
cluster sizes much larger than the critical cluster, but it ignores effects of 



15 



the surface tension and of the local environment (the diffusion field). Dur- 
ing the phase transformation of the (Ising) square lattice-gas model in the 
multi-droplet regime the clusters are not perfectly circular [9,23], but one can 



define an effective radius r = \Jn/ix. We compute at t\ the microscopic in- 
terface velocity for every supercritical cluster i which has been born at to,i, 
as Vi = (ri(ti) — ri(toj))/(ti — to,i), while excluding very recently nucleated 
clusters (to,* < 0.75ti). Here, t\ is the first time the coverage in the actual 
simulation runs reaches 9 = 0.1, thus minimizing the effect of coalescence (see 
below). The resulting distribution of microscopic interface velocities for each 
set of simulation parameters is roughly bell-shaped and yields the mean and 
standard deviations given in Table 7. The mean microscopic interface velocities 
increase somewhat with increasing supersaturation, and in quite a pronounced 
manner with the diffusion rate. Where data are available, the rates from Seki- 
moto's theory are in good agreement. The standard deviations are typically 
about half of the mean. For the simulations in the absence of diffusion, the 
interface velocity compares quite well with the results from a dynamic mean- 
field approximation for a solid-on-solid interface under the same microscopic 
dynamics [42]. This is also shown in Table 7. 

4-2.2 Nucleation and coalescence rates 

In classical nucleation theory, for a sufficiently low supersaturation and for 
a large critical nucleus, nucleation kinetics can be described by the Fokker- 
Planck equation, as the superposition of a drift and a diffusional motion in 
particle-size space [43]. The former is driven by the derivative of the cluster 
free energy F(n) with respect to the cluster size n, and the latter is a stochastic 
fluctuation of the particle size. 

A nucleation event happens when a cluster grows to a size larger than a 
certain cutoff n c . Likewise, a supercritical cluster dissolves when it shrinks 
back to a size less than or equal to n c . When one supercritical cluster splits 
into two supercritical clusters, this is not counted as nucleation, and when 
two supercritical clusters coalesce, this is not counted as dissolution of the 
vanishing cluster. At the top of the free-energy barrier, when n c = n*, these 
processes are dominated by the diffusional motion in n-space. Far away from 
the barrier, this diffusional motion can be neglected. The diffusion length in 
n-space is [43] 



where D = k^T and F(n) is the free energy of a cluster of size n. From 
a parabolic fit of the barrier region of F(n) for u — /i = 0.4 this yields 
l n = 11. In Fig. 10 we compare the raw nucleation and dissolution rates for 
n c = n* = 18 and for n c = n* + l n = 29. The difference between the raw 





(20) 
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nucleation and dissolution rates yields the net nucleation rate. Choosing the 
larger cutoff, the diffusional motion is somewhat reduced, but not eliminated, 
and the net nucleation rate takes longer to reach its plateau value. We thus use 
n c = n* in the following. We further note that the conditions for the validity of 
a Fokker-Planck equation are not strictly fulfilled with the model parameters 
used. 

The nucleation rates obtained with this method are real nucleation rates, I r . 
Figure 11 shows the net nucleation rate for fi — /x = 1.2 in the absence of 
diffusion and for fast diffusion rates. In both cases the net real nucleation rate 
rises steeply in the initial phase and goes through a maximum; it then de- 
creases continuously to zero. We show this for all diffusion rates investigated 
at this supersaturation in Fig. 12, where the time scale is normalized to the 
metastable lifetime, separating changes in the dynamics from a mere acceler- 
ation of the phase transformation. The decrease is much faster when diffusion 
rates are high. This gives an a posteriori justification for using a decaying 
extended nucleation rate in the analysis of the relaxation function at least on 
a qualitative level. However, it is quite noticeable that even in the absence 
of diffusion the real nucleation rate reaches zero when only about half of the 
surface is covered. This is clearly in contradiction to the assumed constant 
extended nucleation rate in KJMA theory (which would result in a drop of 
the real nucleation rate to half its maximum value at this coverage). 

Figure 11 also shows the net coalescence rate, which is the difference between 
the rates at which supercritical clusters vanish by coalescence and appear by 
dissociation. For fast diffusion rates it increases more strongly in the initial 
phase, but saturates earlier and at lower values than in the absence of diffu- 
sion. We show this for all diffusion rates investigated at this supersaturation 
in Fig. 13, where the time scale is normalized to the metastable lifetime. We 
believe that this pattern of the coalescence rate for fast diffusion is related to 
the changes in morphology (Section 4.1.2). The initially more effective coales- 
cence results in a reduction of the density of mid-sized clusters and an increase 
of the correlation length, such that later on fewer clusters are prone to coalesce 
in the more strongly ordered adsorbate layer. The difference between the net 
nucleation and coalescence rates is the net rate at which supercritical clusters 
are formed. Its integral over time perfectly reproduces the time-dependent 
densities of supercritical clusters, corroborating the accuracy of our cluster- 
counting method and analysis. 

In the absence of diffusion, the decrease of the real nucleation rate comes from 
the reduction of the free surface area. Its maximum value at low coverage thus 
gives an estimate of the extended nucleation rate / in the KJMA theory. The 
results are shown in Table 6. They are of the same order as the results from 
Sekimoto's theory, though systematically a little higher. 
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When the microscopic net nucleation rate reaches its maximum, the net coa- 
lescence rates are still close to zero. This implies that the real nucleation rates 
at early times so obtained should be a good estimate for the extended nucle- 
ation rate, since the microscopic nucleation rate is reasonably well separated 
from coalescence processes. Also, the microscopic nucleation rate should be 
very close to the time derivative of the density of supercritical clusters in the 
initial period of the phase transformation. This is indeed borne out by our 
results which are shown in Table 6. 



4.2.2.1 Simulations of the metastable phase Further insight into the 
nature of the nucleation process is gained by simulations in which we suppress 
the decay of the metastable phase and the adsorption of supercritical clusters. 
We start the simulations with an empty surface and at time t = we switch on 
a supersaturation favorable for adsorption. We follow the sizes of all clusters 
at every Monte Carlo step. Whenever a step produces a supercritical cluster, 
this whole cluster is desorbed and the event counted as a nucleation event. 
The histogram for the time between two nucleation events, t nuc \, is shown in 
Fig. 14 as a log-lin plot for fj, — fj, — 1.2; the plots for the other supersatura- 
tions look similar. The probability distribution of i nuc i decays exponentially, 
as expected for a Poisson process. A further indication is the fact that the 
standard deviation of t nuc i equals its mean (£ nuc i) (see Table 6). It must be 
noted that, since there are never any stable clusters present on the surface the 
effects of spatial correlation on the nucleation rate are neglected. The nucle- 
ation rates obtained from the reciprocal of L 2 (i nuc i) are about a factor two 
larger than the microscopic rates obtained from the simulations of the phase 
transformation. This corresponds to a transmission factor of about 1/2 for the 
diffusional motion across the nucleation barrier in n-space. For fast diffusion 
rates the discrepancy is a little larger, which might be an effect of the absence 
of exclusion zones in which no nucleation can take place in these "simulations 
of the metastable phase" . We emphasize that no real metastable phase is pro- 
duced by this method since the complete desorption of supercritical clusters 
disturbs the quasi-equilibrium. 

4-3 Accuracy of the continuum description 

The time-dependent law for the extended nucleation rate that we used for 
the analysis of the phase-transformation dynamics is a mean-field approxima- 
tion whose accuracy is difficult to predict (see Section 2.1). Our simulations 
have shown a strong decrease of the real nucleation rates in the presence of 
fast diffusion. However, not even in the absence of diffusion does the real 
nucleation rate follow exactly the assumptions of the theory (Section 4.2.2). 
Moreover, the two-point correlation functions of the adsorbate phase indicate 
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that the condition of spatially random nucleation is violated when diffusion 
rates become noticeable (Section 4.1.2). A check of the validity of Avrami's 
law applying an exponentially decaying extended nucleation rate to describe 
the phase transformation in the present microscopic model, is to compare the 
kinetic parameter /Wmodei from the fit of the relaxation function to the theory 
with /^micro calculated from the microscopic rates / and v. We assume here 
that the latter represents the correct value. This is shown in Fig. 15. 



In the absence of diffusion and for low diffusion rates the agreement is very 
good. When the diffusion rate increases, we find systematic deviations: the 
result from the fit is increasingly too low. Still the fits yield the correct order 
of magnitude for Iv 2 . Hence, the agreement between the KJMA theory for 
exponentially decaying nucleation rates and the microscopic results is quite 
satisfactory, keeping in mind that the former is only a coarse-grained con- 
tinuum approximation of our microscopic model system. In particular, the 
incomplete separation of the time scales for the formation and the decay of 
the metastable phase is a potential source for deviations and limits the qual- 
ity of the curve fitting. Introducing an additional fitting parameter that shifts 
the time scale might improve the ability of the model to fit the shape of the 
phase-transformation curves, but its physical meaning would be questionable. 



In the limit of instantaneous nucleation, Iq/X should equal the density of nu- 
clei. This can be checked using the microscopic I T at early times as an estimate 
for Iq and the maximum density of supercritical clusters as an estimate for 
the density of nuclei. The results are presented in Table 8. They show quite 
good agreement and the correct trend with the supersaturation, despite the 
deviations of Iv 2 under the respective conditions. Hence, A is indeed physically 
meaningful and not a mere model parameter. 



We point out that the presented theory does not account for spatial corre- 
lations between the nuclei that are a consequence of the diffusion-induced 
exclusion zones. It is difficult to anticipate the consequences of an appropri- 
ate extension of the theory. However, since these correlations shift the overlap 
between growing extended clusters to later times, accelerating the phase trans- 
formation at early times, one might speculate that the discrepancies in Iv 2 
would rather increase. We conclude that even though some of the assumptions 
behind the KJMA theory are not strictly fulfilled, it still provides a reason- 
able description of the phase-transformation kinetics of the Ising lattice gas. 
The effects of adsorbate diffusion can in part be accounted for by assuming a 
decaying nucleation rate. This effectively yields a crossover from continuous 
to instantaneous nucleation when diffusion is the predominant process. 
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4-4 Onset of the crossover 



Our previous results have shown that a crossover of the nucleation from con- 
tinuous to instantaneous sets in when the fit of the relaxation function to 
Eqs. (4), (6), and (15) results in a parameter At of the order of one. Then, we 
find changes in the morphology of the adsorbate phase and a faster decrease 
of the microscopic nucleation rate. As commonly accepted in vacuum sur- 
face science, diffusion of adsorbed particles results in exclusion zones around 
existing clusters, in which nucleation is suppressed, since free particles rather 
aggregate to the cluster than form a new nucleus. Our hypothesis is that these 
exclusion zones become space-filling at early times for high diffusion rates [9]. 
Here we present a rough estimate of when the crossover sets in. 

The average distance between independent droplets is roughly [23] 

R « vt. (21) 

From the competition between desorption and diffusion, a monomer after the 
time it adsorbs performs on average 

jV(oo) = ^ (22) 

Pdes 

steps. Here, from Eq. (10) p&f = ^oexp(— A^/Zcb^) is the diffusion probabil- 
ity of a free monomer, and pd cs = ^exp(— A^/des/feT 1 ) exp(— fx/2 k^T) with 
fx = (/i — /xq) — 2$ is its desorption probability. In a two-dimensional random 
walk the particle travels a distance ^Jn(oo). We compare yOv '(oo) and Rq 
in Table 9. For Adif = 15 the monomer performs on average less than one 
step. The crossover to instantaneous nucleation is complete when the diffu- 
sion length of a monomer is about half the average cluster distance. Despite 
the crudeness of the approximation it gives a satisfactory explanation for the 
onset of the crossover. It correctly reproduces the trends with supersaturation 
and diffusion barrier. 



5 Summary and Conclusions 

We have performed kinetic Monte Carlo simulations of the adsorption dynam- 
ics of a monolayer in a square lattice-gas model using Arrhenius dynamics with 
local barriers. The model Hamiltonian is grand-canonical, meaning that local 
adsorption/desorption equilibrium is approximately established. This makes 
it suitable for submonolayer electrochemical deposition under control of the 
potential difference between solution and electrode and when mass-transport 
limitations to the surface are insignificant. The elementary steps included in 
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the model are adsorption, desorption, and diffusion; our main interest is in 
the influence of the latter on the phase-transformation dynamics. The dif- 
fusion rates can be controlled by the barrier height for the diffusion steps. 
We have obtained results for the phase-transformation kinetics, the adsorbate 
morphology, and the microscopic rates of nucleation, cluster growth (interface 
velocity), and coalescence. For our model calculations we have used an ex- 
tension of Bortz et al.'s n-fold way algorithm that contains both adsorption 
(spin-flip) and diffusion (spin-exchange), and whose details we present in this 
paper. 

The nucleation rate and the interface velocity and consequently the phase- 
transformation kinetics are accelerated by both an increase of the supersatu- 
ration (the driving force) and the diffusion rate. The higher the diffusion rate, 
the more likely it is that a walker aggregates with sub- or supercritical clusters 
to participate in nucleation or growth, respectively. 

Around supercritical clusters, diffusion causes a local reduction of the nucle- 
ation rate in a zone with the extension of the monomer diffusion length, since 
it is more likely that monomers reach the existing cluster than that they cre- 
ate a new one. This causes a reduction of the cluster density and an increase 
of the correlation length. Moreover, diffusion changes the spatial correlation 
between the nuclei, which is no longer random, and creates minima in the 
two-point correlation functions. 

For a theoretical description of the phase transformation we use the Kolmogorov- 
Johnson-Mehl-Avrami (KJMA) theory as a coarse-grained approximation to 
our microscopic model system. The theory is based on the assumption that 
nucleation is spatially random. Using an appropriate law for the time depen- 
dence of the nucleation rate, one can obtain an estimate for Iv 2 from a fit 
of the time dependence of the relaxation function of the metastable phase, 
but no separate estimates for the nucleation rate / and the interface veloc- 
ity v. In the absence of diffusion, this gives very accurate results assuming a 
constant extended nucleation rate / (continuous nucleation), even though our 
microscopic results suggest that the real nucleation rate decreases somewhat 
too fast with time for this assumption to be strictly fulfilled. With diffusion, 
though, the exclusion zones around supercritical clusters grow with the cluster 
radius, such that the real nucleation rate diminishes more strongly with time 
than expected from the decrease of the free surface fraction (and than in the 
case without diffusion). Assuming an extended nucleation rate that decreases 
exponentially with time from an initial value I , we can extend the analysis to 
this case and still extract an estimate for IqV 2 . This estimate is less accurate, 
but still correct by order of magnitude. We emphasize that correlation effects, 
which are also important in these cases, are neglected, since we use Avrami's 
law. When diffusion is predominant, the exclusion zones become space-filling 
at early times of the phase transformation, and the dynamics resemble in- 
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stantaneous nucleation. A fit to the respective form of Avrami's law would, 
however, not yield an estimate for IqV 2 . Instead, IqV 2 /X is obtained from this 
analysis, where A is the time constant of the decay of /. Iq/\ is the density of 
nuclei. We emphasize that nucleation is still homogeneous. Hence, in the elec- 
trodeposition of a two-dimensional film that proceeds by a first-order phase 
transition, instantaneous nucleation need not imply heterogeneous nucleation 
at defects. Also, intermediate cases are possible that would lead to non-linear 
Avrami plots and non-integer Avrami exponents. A strong indication for our 
model to be a good description of a given experimental system could be found 
in the two-point correlation functions, which develop minima for an Avrami 
exponent less then 3. They can be obtained from scattering experiments [44] 
or by analysis of microscopic images. 

In the absence of diffusion, it is possible to obtain separate estimates for I 
and v from the variance of the coverage vs. time using Sekimoto's theory for 
the two-point correlation function. This theory is based on spatially random 
nucleation as well, and the analysis breaks down when diffusion rates become 
noticeable and a minimum in the correlation function appears. 

The time a walker spends on the surface before it desorbs increases with 
higher supersaturation, such that in this case the influence of diffusion on the 
dynamics of the phase transformation at the same rate is stronger. Hence, the 
crossover to instantaneous nucleation is more complete for the same diffusion 
barrier when the supersaturation is greater. 

In a previous work [9] we have presented kinetic Monte Carlo simulations of 
the same model using Glauber dynamics, but without varying the supersat- 
uration, and performing a less stringent analysis of the data. Still, we found 
from the time- dependence of the relaxation function qualitatively the same 
crossover to instantaneous nucleation as in the present investigation. In vari- 
ous studies for solid-on-solid (SOS) and Ising models, a comparison of different 
transition probabilities in kinetic Monte Carlo dynamics has shown that the 
stochastic dynamics can have a strong influence on the microscopic properties 
of the evolving system, namely interfacial structure, nucleation rates, and in- 
terface velocities (see, e.g., Refs. [42,45,46]). Distinctions are made between 
so-called "soft" dynamics, in which the field (supersaturation) and coupling 
(interaction) constant terms in the transition probability factorize, and "hard" 
dynamics, in which they do not. The dynamics in the present work are soft, 
while in Ref. [9] they are hard. Another difference is that here we use local 
barriers in contrast to Ref. [9]. In both investigations, the stochastic dynamics 
contain conserving and non-conserving steps with respect to the order param- 
eter (the coverage). This paper does not aim at a systematic investigation 
of the influence of the stochastic dynamics; in fact the diffusion is treated in 
such different ways that a quantitative comparison is difficult. It is, however, 
possible to conclude that the decrease in local nucleation rates around clus- 
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ters and the concomitant development of spatial correlations translate into 
qualitatively similar changes in the phase-transformation dynamics. 
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Appendix 

Here we introduce our single-site n-fold way scheme for the square lattice- 
gas model with both adsorption/desorption and diffusion steps. For compu- 
tational convenience, we present the classification scheme in Ising language, 
s = +1 representing a spin up (occupied site) and s = — 1 representing a spin 
down (empty site). For the sake of consistency, we also present the energy 
differences in terms of the Ising Hamiltonian. Within Ising terminology, ad- 
sorption/desorption is equivalent to spin flip and diffusion to spin exchange. 
The former involves only one particle, while the latter involves a pair. Hence, 
care must be taken that diffusion steps are not counted twice when we set up 
our list of all elementary steps in the lattice. For each particle, only spin ex- 
change with two of the four nearest neighbors in orthogonal directions (here: 
right and below, respectively) is to be considered. The spin exchange with the 
nearest neighbors in the opposite directions (left and above, respectively) is 
included in the neighbors' respective set of moves. In lattice-gas terms, instead 
of the diffusion of a particle to the left or upwards, we consider the diffusion 
of a vacancy to the right or downwards, starting from the respective nearest- 
neighbor position. We use the following labels for the three elementary steps 
associated with each lattice point: m = 1 (adsorption/desorption), m = 2 
(diffusion to [spin exchange with] the nearest-neighbor position to the right, 
and m = 3 (diffusion to [spin exchange with] the nearest-neighbor position 
below). We call the classes for the adsorption/desorption steps particle classes 
and for the diffusion steps pair classes. For the particle classes, we need to 
consider the value (spin) of the central particle C and of its four nearest 
neighbors, s&. For the pair classes, we need to consider the spins s$ and Sj of 
two neighboring particles, C and C, respectively. C and C form the pair S. 
Moreover, we consider the spins Sk of all nearest neighbors of C, and the spins 
si of all the nearest neighbors of C. The central particle C is the one to which 
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the respective pair class belongs. C is called the exchange partner. 



For the particle classes c a /d, the scheme is shown in Table A-1 [32]. From it 
follows the simple relation 

c a /d = ^ - 5si - s kj • (A" 1 ) 

The difference of the configurational energy resulting from the respective ad- 
sorption/desorption step is obtained from 

H(b) -H(a) = AH = 2J Si J2sk + 2H Si . (A-2) 

k 

For the pair classes c^a, only the cases where Sj ^ Sj need to be considered; 
the scheme is shown in Table A-2. It follows the simple relation 

Cdif = -ySi + 2 X> fc + ^I>i + y. (A-3) 

The difference of the configurational energy performing the respective diffusion 
step is obtained from 



H(b) - H{a) = AH = J 



V k i J 



(A-4) 



Not only are the classes changed by moves in which the respective site is 
involved, but also by moves in the neighborhood. The classes of up to some of 
the 4th-nearest neighbors have to be updated after a move. See Fig. A-1 for 
an overview of the sites affected by a move, and Fig. A-2 for the designation 
of the neighboring positions. We denote the spin of an affected site with U 
when it is the central particle X, and with tj when it is the exchange partner 
X' of X. X and X' form the pair T. The spins of all the nearest neighbors of X 
and X' are labeled tk and ti, respectively. Tables A-3, A-4 and A-5 summarize 
how the classes in the neighborhood are affected. 
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Figure captions 

Fig. 1. Fraction of diffusion steps as a function of the coverage during adsorption of 
a monolayer. Diffusion barrier A d if = 15. 

Fig. 2. Acceleration of the phase transformation with increasing electrochemical 
potential. Adif = oo. 

Fig. 3. Acceleration of the phase transformation with diffusion, jjl — liq = 0.4. 

Fig. 4. Crossover from continuous to instantaneous nucleation for li — ho = 1.2. (a) 
Logarithmic plot of the relaxation function vs. the cube of the normalized time, as 
expected for continuous nucleation, Eqs. (2) and (4). (b) Same vs. the square of 
the normalized time, as expected for instantaneous nucleation, Eqs. (3) and (4). In 
both cases, the straight line is a guide to the eye. 

Fig. 5. Fits of the relaxation function to Avrami's law for exponentially decaying 
nucleation rate. Dotted lines: simulations. Continuous lines: fits. The curves are, 
from right to left, for A^f = oo, 15, 10, 5. \i — ho = 1.2. 

Fig. 6. Two-point correlation functions G(r) as functions of the distance r for various 
coverages, (a) ll~ llq = 0.4, A di f = oo. (b) // — Ho = 1.6, A^f = 5. Note the different 
scale along the horizontal axis in the two parts. 

Fig. 7. Correlation length / vs. 8, estimated from the inverse of the initial slope 
of the normalized correlation function G(r). Thin lines: fj, — /j,q = 0.4. Bold lines: 
M - Mo = 1-6. 

Fig. 8. Cluster-size distribution for jjl — fj-o = 1.6 at coverage 8 = 0.3. The vertical 
line marks the critical cluster size n* . 

Fig. 9. Snapshots of the adsorbate phase during the phase transformation. From left 
to right: fj, - fi = 0.4, A dif = oo; fx - fj, = 0.4, A dif = 5; li - li = 1.6, A dif = oo; 
li — no = 1.6, Adif = 5. Upper row: 8 = 0.3. Lower row: 8 = 0.7. 

Fig. 10. Nucleation and dissolution rates and net nucleation rate for a cutoff size 
equal to the critical cluster size (thin lines) and one diffusion length in n-space 
larger (bold lines), li — Ho = 0.4, A d if = 15. 
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Fig. 11. Net nucleation and coalescence rates and net cluster formation rate, (a) 
Adif = oo. (b) Adif = 5. In both cases, jj, — /j,q = 1.2. Note the different scales in the 
two parts. 



Fig. 12. Net nucleation rate for various diffusion barriers, fi — fio = 1.2. 



Fig. 13. Net coalescence rate for various diffusion barriers, n — no = 1.2. 



Fig. 14. Probability distribution of the time between two nucleation events in the 
simulations of the metastable phase. /U — /io = 1.2. 



Fig. 15. The kinetic parameter Iv 2 from the microscopic rates / and v and from the 
fit of the relaxation function to the theory. 



Fig. A-l. (Color online) An overview of the affected classes of sites by a move of the 
central particle, (a) particle classes affected by adsorption/desorption: central site 
(red) and nearest neighbors (green), (b) particle classes affected by diffusion: the 
switching pair (red) and its nearest neighbors (green), (c) pair classes affected by 
adsorption/desorption: the pairs that contain the central particle (red) and those 
that contain one of its nearest neighbors (blue), (d) pair classes affected by diffusion: 
the diffusing pair (yellow), pairs that contain the central particle (red), pairs that 
contain the exchange partner (brown), pairs that contain nearest neighbors of the 
central particle (blue), of the exchange partner (pink), and of both (green). Bold 
arrows: the pair starts at the respective particle or nearest-neighbor site; thin arrows: 
the pair points to the respective site. 



Fig. A-2. Designation of the affected nearest (NN), next-nearest (NNN), 3rd-nearest 
(3NN), and 4th-nearest (4NN) neighbor position when the central site undergoes a 
move. The numbers are the labels used in Tables A-3, A-4, and A-5. 
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Tables 



Table 1 

The average fraction of diffusion steps throughout the phase transformation. 

5 



M - Mo 


Adif = 15 


Adif = 10 


Adif = 


0.4 


0.221 


0.817 


0.986 


0.8 


0.239 


0.830 


0.987 


1.2 


0.247 


0.834 


0.987 


1.6 


0.250 


0.834 


0.986 



Table 2 

The critical cluster size n* as a function of the electrochemical potential. 

H - Mo n* 
0.4 18 
0.8 7 
1.2 4 
1.6 3 

Table 3 

The metastable lifetime r in MCSS. r is the time when the relaxation function 
4> = 1/2. 

5 



M - Mo 


Adif = oo 


A di f = 15 


Adif = 10 


Adif = 


0.4 


270000 


260000 


216000 


137000 


0.8 


59000 


56800 


49200 


38800 


1.2 


27500 


26800 


24100 


21300 


1.6 


16500 


16200 


15150 


14300 



Table 4 

The kinetic parameter Iqv 2 in MCSS -3 , from a fit of the relaxation function (f> to 
Avrami's law with an exponentially decaying nucleation rate. 

M - Mo A dif = oo A dif = 15 A dif = 10 A dif = 5 

0.4 3.0 x 10" 17 3.3 x 10" 17 5.9 x 10~ 17 3.3 x 10~ 16 

0.8 2.7 x 10~ 15 3.0 x 10~ 15 5.0 x 10~ 15 1.9 x 10~ 14 

1.2 2.6 x 10" 14 2.8 x 10~ 14 4.5 x 10~ 14 1.2 x 10~ 13 

1.6 1.2 x 10" 13 1.2 x 10~ 13 1.8 x 10~ 13 4.0 x 10~ 13 
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Table 5 

The parameter At, from a fit of the relaxation function cf> to Avrami's law with 
an exponentially decaying nucleation rate. At indicates how much the extended 
nucleation rate I has decayed at the metastable lifetime. Standard errors are given 
in brackets below the values. 



At - Mo 


Adif = oo 


A dif = 15 


Adif = 10 


Adif = 


0.4 


6.7 x 10~ 3 


1.1 x 10" 4 


2 x 10~ 2 


1.70 




(6.6 x 10~ 3 ) 


(1 x 10~ 5 ) 


(1 x 10- 2 ) 


(0.04) 


0.8 


8.3 x 10~ 5 


7 x 10~ 5 


0.36 


3.52 




(6 x 10~ 6 ) 


(1 x 10" 5 ) 


(0.02) 


(0.03) 


1.2 


5.8 x 10~ 5 


3 x 10~ 2 


0.65 


4.13 




(8 x 10~ 6 ) 


(2 x 10~ 2 ) 


(0.01) 


(0.02) 


1.6 


6.9 x 10~ 5 


1.5 x 10~ 4 


0.79 


4.20 




(8 x 10~ 6 ) 


(3 x 10~ 5 ) 


(0.01) 


(0.01) 



Table 6 

Nucleation rate per site and MCSS. Extended nucleation rate / from fit of I? var(#) 
to Sekimoto's theory; real nucleation rate I r from the maximum of the time- 
dependent net nucleation rate obtained from cluster counting; I T from the ini- 
tial slope of the density of supercritical clusters; and I r from simulations of the 
metastable phase. For the simulations of the metastable phase, also the mean and 
standard deviation of the time in MCSS between nucleation events for the lattice 
size are given. For details of the methods, see the text. 

\i — no Adif Sekimoto microscopic cluster density metastable 





I 




Ir 




Ir 




Ir 




(^nucl) 


a 


0.4 oo 


7.2 x 10" 


-9 


1.4 x 10" 


-8 


1.2 x 10" 


-8 


2.8 x 10- 


-8 


556.9 


558.3 


15 


8.3 x 10" 


-9 


1.3 x 10' 


-8 


1.1 x 10" 


-8 


2.9 x 10- 


-8 


528.3 


523.4 


10 


1.1 x 10' 


-8 


1.8 x 10' 


-8 


1.4 x 10' 


-8 


3.8 x 10- 


-8 


402.1 


400.6 


5 






3.0 x 10" 


-8 


2.5 x 10" 


-8 


8.7 x 10" 


-8 


176.1 


174.8 


0.8 oo 


1.6 x 10- 


-7 


2.3 x 10" 


-7 


2.1 x 10" 


-7 


4.6 x 10" 


-7 


33.2 


33.2 


15 


2.0 x 10" 


-7 


2.5 x 10" 


-7 


2.3 x 10" 


-7 


4.9 x 10" 


-7 


31.2 


31.3 


10 






3.1 x 10" 


-7 


2.7 x 10- 


-7 


6.6 x 10" 


-7 


23.1 


23.3 


5 






4.3 x 10" 


-7 


3.9 x 10" 


-7 


1.2 x 10- 


-6 


12.8 


12.6 



continued on next page 
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Table 6 
continued 

fj, — no Adif Sekimoto microscopic cluster density metastable 





/ 




It 


It 


It 


Vermel / 


a 


1.2 oo 


6.4 x 10" 


_7 


9.7 x 10~ 7 


9.1 x 10~ 7 


1.8 x 10~ 6 


8.5 


8.5 


15 


8.3 x 10" 


-7 


1.1 x 10~ 6 


9.5 x 10~ 7 


1.9 x 10~ 6 


8.0 


7.9 


10 






1.2 x 10~ 6 


1.1 x 10~ 6 


2.5 x 10~ 6 


6.1 


6.1 


5 






1.5 x 10~ 6 


1.3 x 10~ 6 


3.9 x 10~ 6 


4.0 


3.9 


1.6 oo 


1.8 x 10" 


-6 


2.3 x 10~ 6 


2.1 x 10~ 6 


3.9 x 10~ 6 


4.0 


3.9 


15 


1.8 x 10" 


-6 


2.4 x 10~ 6 


2.2 x 10~ 6 


4.1 x 10~ 6 


3.7 


3.7 


10 






2.7 x 10~ 6 


2.4 x 10~ 6 


5.1 x 10~ 6 


3.0 


3.0 


5 






2.9 x 10~ 6 


2.6 x 10~ 6 


7.0 x 10~ 6 


2.2 


2.2 



Table 7 

Interface velocity v in unit length per MCSS. From fit of L 2 var(0) to Sekimoto's 
theory; from the increase of the effective radius of the clusters during the initial 
phase of the simulation (standard deviation in parentheses); solid-on-solid interface 
from a dynamic mean-field approximation [42]. For details of the methods, see the 
text. 



At - Mo 


Adif 


Sekimoto 


microscopic 


(standard deviation) 


MFA 


0.4 


oo 


6.5 x 10~ 5 


5.3 x 10- 5 


(2.8 x 10~ 5 ) 


4.26 x 10~ 5 




15 


6.4 x 10~ 5 


5.5 x 10~ 5 


(2.8 x 10~ 5 ) 






10 


7.4 x 10~ 5 


7.0 x 10~ 5 


(3.3 x 10~ 5 ) 






5 




1.4 x 10~ 4 


(5.5 x 10~ 5 ) 




0.8 


oo 


1.3 x 10~ 4 


1.3 x 10~ 4 


(8.0 x 10~ 5 ) 


8.58 x 10~ 5 




15 


1.2 x 10~ 4 


1.3 x 10~ 4 


(8.1 x 10~ 5 ) 






10 




1.8 x 10~ 4 


(9.9 x 10~ 5 ) 






5 




4.0 x 10~ 4 


(1.6 x 10~ 4 ) 




1.2 


oo 


2.0 x 10~ 4 


1.8 x 10~ 4 


(1.3 x 10~ 4 ) 


1.30 x 10~ 4 




15 


1.8 x 10~ 4 


2.0 x 10~ 4 


(1.4 x 10~ 4 ) 






10 




2.9 x 10~ 4 


(1.7 x 10" 4 ) 






5 




6.7 x 10~ 4 (2.8 x 10~ 4 ) 
continued on next page 
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Table 7 
continued 

fi — fio Adif Sekimoto microscopic (standard deviation) MFA 

1.6 oo 2.5 x 1(T 4 2.3 x (1.9 x 10"" 4 ) 1.76 x 10~ 4 

15 2.6 x 10" 4 2.5 x 10" 4 (1.9 x 10" 4 ) 

10 3.9 x 1(T 4 (2.4 x 10~ 4 ) 

5 9.4 x 1(T 4 (3.8 x 10~ 4 ) 



Table 8 

Ratio of nucleation rate Iq and fitting parameter A as estimate for the density of 
nuclei, which is approximated by the maximum density of supercritical clusters. 
Adif = 5. 

fi-fio VA max(p supcrcrit ) 



0.8 


4.7 x 10" 


-3 


3.1 x 10" 


-3 


1.2 


7.5 x 10" 


-3 


5.4 x 10' 


-3 


1.6 


9.7 x 10" 


-3 


7.1 x 10" 


-3 



Table 9 

Estimate of the onset of the crossover in the phase-transformation dynamics. Dif- 
fusion steps of a monomer before desorption A r (oo); diffusion length of a monomer 
Id = a/^V(oo); average distance between clusters Rq. 

H~Ho A dif N(oo) l d Ro k/Ro 

0.4 15 14.3 

10 1.9 1.4 15.1 0.09 

5 30.4 5.5 19.2 0.29 

0.8 15 7.4 

10 2.2 1.5 8.9 0.17 

5 34.0 5.8 15.5 0.37 

1.2 15 5.4 

10 2.4 1.6 7.0 0.23 

5 38.0 6.2 14.3 0.43 

1.6 15 4.1 

10 2.7 1.6 5.9 0.27 

5 42.4 6.5 13.4 0.49 
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Table A-l 

The particle classes for the adsorption/desorption step 



s i J2k S k 


c a/d 


AH 


1 4 


1 


2H + 8J 


2 


2 


2H + AJ 





3 


2H 


-2 


4 


2H -AJ 


-4 


5 


2H - 8 J 


-1 4 


6 


-2H-8J 


2 


7 


-2H - 4J 





8 


-2H 


-2 


9 


-2H + 4J 


-4 


10 


-2H + 8 J 



Table A- 2 

The pair classes for the diffusion step 

s i s 3 J2k S k J2l s l c dif AH 







-2 


1 








2 


-4J 


2 


3 


-8J 


4 


4 


-12J 


-2 


5 


4J 





6 





2 


7 


-4J 


4 


8 


-8J 


-2 


9 


8J 





10 


4J 


2 


11 





4 


12 


-4J 


-2 


13 


12J 





14 


8J 


2 


15 


4J 


4 


16 






continued on next page 
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Table A-2 
continued 





J2iSi 


Cdif 


AH 


-2 


-4 


17 







-2 


18 


4J 







19 


8J 




2 


20 


12J 





-4 


21 


-4J 




-2 


22 










23 


4J 




2 


24 


8J 


2 


-4 


25 


-8J 




-2 


26 


-4J 







27 







2 


28 


4J 


4 


-4 


29 


-12J 




-2 


30 


-8J 







31 


-4J 




2 


32 
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Table A-3 

Update of the particle classes c a /d after a move: c ncw = c id + Ac. For the designation 
of the spins and particles, see the text. For the meaning of the symbols for the site, 
see Fig. A-2. 

move description Ac site 
a/d X is C +5sj j0 id central 

X is NN of C + Si oW NN 0, 1, 2, 3 



dif X is C +4si j0 id central * t 

X is C +4s ii0 id NN 1 * 

NN 3 1 

NN of X is C +s ii0 id NN 0, 2, 3 * 

NN 0, 1, 2 t 
NN of X is C +s ii0 id NNN 1, 2; 3NN 1 ; 



NNN 1, 3; 3NN 3 f 



* move was diffusion to the right: C is NN 1 
t move was diffusion to site below: C is NN 3 



Table A-4 

Update of the pair classes after an adsorption/desorption move: c new = c Q id + Ac. 
m = 2: affected class is Cdif for diffusion to nearest-neighbor position to the right 
(NN 1); m = 3: affected class is Cdif for diffusion to nearest-neighbor position below 
(NN 3). For the designation of the spins and particles, see the text. For the meaning 
of the symbols for the site, see Fig. A-2. 

site 

description Ac 

m = 2 m = 3 

X is C — 32sj i0 idtj central central 

X'isC -32*^-16*^ NN NN 2 

NN of X is C -4a i)OM NN 1, 2, 3 NN 0, 1, 3 

NN of X' is C -s ii0 i d NNN 0, 3; 3NN NNN 0, 2; 3NN 2 
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Table A-5 

Update of the pair classes after a diffusion move: c new = c Q id + Ac. m = 2: affected 
class is Cdif for diffusion to nearest-neighbor position to the right (NN 1); m = 3: 
affected class is Cdif for diffusion to nearest-neighbor position below (NN 3). For the 
designation of the spins, particles and pairs, see the text. For the meaning of the 
symbols for the site, see Fig. A-2. 

site 

description Ac 







m = 2 


m = 3 


T is S 


+16sj i0ki 


central * 


central + 


X is C, T _L S 


— 32sj j0 idtj + 4s ij0l d 


central ' 


central * 


X' is C 


— 32sj j0 idtj — 15sj i0 id 


NN0*t 


NN 2* t 


X is C 


-S2sj :0 idtj + 4sj 


NN 1* 


NN 1* 






NN 3 1 


NN 3 1 


X' is C 


~ 32sj i0 id£i — 15sj j0 id 


NNN 3 1 


NNN 2* 


T || S 


— 3sj id 


NN 2, 3* 


NN 0, 1 1 


NN of X is C 


-4si, id 


NN 1,2 + 


NN 0, 3 * 


NN of X' is C 


Sj j0 id 


NNN 0, 3; 3NN * 


NNN 1; 3NN 1* 






NNN 0; 3NN + 


NNN 1, 3; 3NN 3 1 


NN of X is C 


Sj j0 ld 


NNN 1, 2; 3NN 1* 


NNN 0; 3NN 2 * 






NNN 1; 3NN 3 1 " 


NNN 0, 2; 3NN 2 1 


NN of X' is C 


Sj j0 ld 


4NN 0, 3 1 


4NN 1, 2* 



* move was diffusion to the right: C is NN 1 
t move was diffusion to site below: C is NN 3 
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Fig. 9. 
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Fig. A-l. (greyscale) 
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Fig. A-l. (color) 
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Fig. A-2. (greyscale) 
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Fig. A-2. (color) 
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